Relativistic Astrophysics with Resonant Multiple Inspirals 
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We show that a massive black hole binary might resonantly trap a small third body {e.g. a neutron 
star) down to a stage when the binary becomes relativistic due to its orbital decay by gravitational 
radiation. The final fate of the third body would be quite interesting for relativistic astrophysics. 
For example, the parent binary could expel the third body with a velocity more than 10 % of the 
speed of light. We also discuss the implications of this three-body system for direct gravitational 
wave observation. 



I. INTRODUCTION 

The orbital period of Pluto is 3/2 times that of Nep- 
tune, and their mutual stability is maintained by this 
simple commensurability (termed the mean motion res- 
onance). In addition to this well-known fact, the Solar 
system has various forms of the mean motion resonances 
[l|. Furthermore, more than 8 exoplanet systems are 
known to have two planets in mean motion resonances. 
One of the notable properties here is that once two plan- 
ets are trapped in a stable resonance relation, they often 
keep the state for a long time. For example, recent nu- 
merical studies showed that some resonant trappings are 
strong enough to be preserved against planet migration 
during which orbital radii of twq^ planets decreased by 
more than 1 order of magnitude [2| . 

Black hole (BH) binaries are considered as fascinating 
sources of broad astrophysical phenomena, and, at the 
same time, their inspirals and mergers are promising tar- 
gets for gravitational wave (GW) observation that would 
also provide us with ideal opportunities to test funda- 
mental aspects of gravity. As in the case of planetary 
systems, a massive BH binary might resonantly trap a 
small body, e.g. through interaction with its surrounding 
disk, and shrink its orbit by dissipative processes. In this 
paper we call this kind of three-body system a resonant 
multiple inspiral (RMI), and study the evolution and the 
final fate of an RMI, including the post-Newtonian ef- 
fects. We find that a BH binary has potential to keep 
a third body down to a stage when the binary becomes 
relativistic due to its orbital decay by gravitational radi- 
ation. We also mention impacts of RMIs on relativistic 
astrophysics and future GW observation. 

II. STABLE EQUILIBRIUIVI POINTS 

We study the evolution of RMI using a circular re- 
stricted three-body problem, namely, analyze the motion 
of a test particle trapped by a BH binary in a circular or- 
bit with masses (1 — IJ,)M and fiM [M: total mass of the 
binary, ^(< 0.5): the mass ratio]. We use the geometri- 
cal unit G ~ c ^ M = 1 with which the Schwarzschild 
radii of the two BHs are 2(1 — /i) and 2/i respectively. If 



necessary, we intentionally recover the mass parameter 
M to show the actual scales of physical quantities. We 
introduce a parameter ri, for the binary separation. 

In this paper, we concentrate on 1:1 mean motion res- 
onances as tractable but intriguing examples (see e.g. Q 
for recent related studies). In this case the test particle 
moves around the equilibrium points L4 or L5. The two 
points are at almost equilateral positions relative to the 
parent binary on its orbital plane. After the pioneering 
work by Lagrange in 1772, various properties have been 
theoretically investigated for the two points. In the real 
Universe, the Sun- Jupiter system {^ ^ lO"'^) has a large 
number of asteroids known as Trojans around its L4 and 
L5 (first discovered in 1906) [lj,i^], and their origin is still 
on active debates [1, Q ■ Similar objects have been found 
e.g. for the Sun-Mars (/z ~ 2 x 10""^) or Saturn- Tethys 
(/i ^ 10~^) systems [l|- Here Tethys is the fifth largest 
moon of Saturn. 

In order to describe the motion of a test particle around 
L4 or L5, it is convenient to introduce a normalized frame 
(Xpf, Y)v, Z]^) that is co-rotating with the parent binary 
around its barycenter. In this frame, the larger BH is at 
(/i, 0,0), the smaller one is at (—(1 — /i),0,0) and their 
separation is unity. The Zn axis is the rotation axis of 
the parent binary and normal to its orbital plane [24]. 
Meanwhile, the positions of the two equilibrium points 
L4 and L5 are given by {Xl, — Yl, 0) and {Xl,Yl,0) with 

x, = -l + , + Sdi|±i^, „(.,.), a, 

where the terms ex r^j" represent the first-order post- 
Newtonian (IPN) corrections [a, Hi- 

We briefiy summarize the results of the linear sta- 
bility analysis for a test particle around the equilib- 
rium points L4 and ^5 [ll, [Sj. In Newtonian mechan- 
ics, small perturbations around these points are given 
by the superposition of two stable oscillating modes for 
^ < ^1 = 1/2 - 769/18 = 0.038521 [|. They are known 
as the epicyclic motion with the frequency uje and the 
libration motion with ujl{< loe) [1[- The two frequencies 



are given by 



^E.L — ^BN 



1 ± y/l - 27m(1 - m) 



(3) 



with the orbital frequency of the parent binary lobn = 
(Af/rj^)^/^ defined at the Newtonian order (correspon- 
dence of signs, uje : + and ujl '■ —) [ll- The two fre- 
quencies degenerate at /i = /ii. For a larger mass ratio 
/i > /ii, the perturbation becomes unstable. 

With IPN analysis, the two basic frequencies ue and 
ojl have the correction terms of 0{r'^ ). It is worth 
stressing that the ratio we/'jJl now depends not only on 
the mass ratio /i but also on the separation rb of the par- 
ent binary. The critical mass ratio ni is given explicitly 
by ^ll = 0.03852 - 0.29056Af/r6 [1]. 



III. EVOLUTION OF RMI 

An interesting question here is how the resonant trap- 
ping changes with the evolution of the parent binary. 
Below, we numerically study this issue, as a restricted 
circular three-body problem. 

For the evolution of the circular parent BH binary, we 
can follow its IPN orbit almost analytically, including the 
orbital decay by the radiation of GW (see the Appendix 
for formulations and basic numerical scheme). Gravi- 
tational radiation reaction generates a dissipative force 
starting at 2.5PN order 9]. Below this order, the sys- 
tem is conservative. Since we assume that the mass of 
the third body (test particle) is negligible, the dissipa- 
tive evolution is controlled by the motion of the parent 
binary, through the fifth time derivative of its quadrupole 
moment. After some algebra, the decrease of the orbital 
separation of the parent binary is described by the stan- 
dard expression ]^ 



dn 
~dt 
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(4) 



Accordingly, the rotation cycle N of the binary before 
the coalescence is estimated by 



7V = 



1 



647r//(l - ^) 






(5) 



Given the positions of the parent binary, we can numer- 
ically follow the evolution of the test particle. Note that 
the dissipative acceleration directly works also on the test 
particle, as easily understood from the fact that the en- 
ergy loss rate is proportional to the square of the GW 
amplitude 9] . This effect might induce interesting effects 
for a 1:1 resonance, due to the apparent phase coherence 
of the three particles. 

Since we cannot properly handle the strong gravity 
around the two BHs with our IPN equation of motion, 
we conservatively terminate our integration when (i) the 
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FIG. 1: Evolution of the orbit of a test particle around the L5 
point in the normalized corotating frame (Xjv, Yn). The mass 
ratio of the parent binary is /i = 0.027. We show the orbits at 
three different separations rb — 200, 95 and 54. Each figure is 
given for ~ 8 rotation cycles of the parent binary. The initial 
conditions are q^ — 0.002, q^ — 0.001 at rj, — 200. The upper 
and bottom figures are shifted toward the vertical direction 
by ±0.004. 



separation of the parent BH binary reaches the inner- 
most stable circular orbit at r;, = 6 [9| or (ii) similarly, 
the test particle goes into 3 Schwarzschild radii of each 
parent BH. As for the initial position and velocity of the 
test particle, we simply introduce two small parameters 
Qx and qz related to its position in the corotating frame 
(Xjv, Yn, Zn) as (Xl(1 + q^), ±YL,2XLqz) with Xl and 
Yl given in eqs.(IT|) and (HJ at the IPN order. Here the 
parameter q^ sets the perturbation of the third body 
within the orbital plane of the parent binary, and an- 
other one qz is related to the inclination with respect to 
the plane. The particle is released with an initial veloc- 
ity with which it is at rest in the corotating frame. Note 
that we made perturbative treatment of the IPN correc- 
tion terms. Thus, even with g^; = g^ = 0, a small initial 
perturbation of 0{r^^) is induced around the equilib- 
rium points due to the truncation effect (see e.g. [101] for 
a similar case). 

We hereafter discuss the results only from the L^ re- 
gion. We have compared the evolution of the test parti- 
cles starting from both the L4 and L5 regions, but have 
not found notable qualitative differences. In Fig.l we 
show the orbits of a test particle at three different epochs. 
This is a typical example for the evolutionary behavior of 
the perturbation around the L5 point. Each figure clearly 
shows the characteristic profile of a tadpole orbit that is 
given by the superposition of the two elliptical motions; 
the short-period epicyclic mode (with the frequency ue) 
and the long-period libration mode (with ojl) [1|. 
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FIG. 2: The unstable separation ri,„ of a parent binary as 
a function of its mass ratio fi. The points are given from 
the evolution of the test particle initially placed around the 
1/5 point at rj — 100. The solid curve is derived from the 
stability condition lul ~ ijJe at IPN and has the asymptotic 
profile rju — > co at /i = 0.038521. The dashed and dot dashed 
curves are given for lje = 2ajL and cue ~ Sljl with the critical 
mass ratios fi = 0.024293 and fi = 0.013516. respectively. 



strong branch around /i ~ 0.02. We found that this 
is caused by a resonant instability due to the coupling 
of the epicyclic and Hbration modes around uje = 2ujl, 
corresponding to the specific mass ratio /i = 0.024293 for 
the Newtonian analysis J13l|. With the PN effects, the ra- 
tio loe/i-^l depends not only on the mass ratio /z but also 
on the binary separation rh, as mentioned before. There- 
fore, even if a parent binary has a mass ratio satisfying 
loeI'^l > 2 at the Newtonian limit r^ — > oo, it could 
match the unstable condition ujeI'-^l = 2 at a finite rf, 
due to the orbital shrink by gravitational radiation. As 
a result, with the general relativistic effects, the binary 
is affected by the resonant instability for a broader mass 
range /i, in contrast with a purely Newtonian system. 

Figure 2 also shows a branch associated with uje = 3wl 
corresponding to the specific mass ratio [i = 0.013516 
for the Newtonian limit [13]. However, this higher-order 
effect is relatively weak, and many binaries can safely go 
through the unstable separation, as shown in Fig. 2. 



IV. FINAL FATE OF RMI 



With a simplified Hamiltonian and an associated adi- 
abatic invariant [ll|, Fleming and Hamilton [l2[ analyt- 
ically predicted how migration of Jupiter affects the or- 
bital evolution of its Trojan asteroids (with Newtonian 
mechanics). Since the time scale of gravitational radia- 
tion is much larger than the orbital period of the parent 
BH binary (at least for r\, given in Fig.l), their analytical 
prediction might fit with our study. Indeed, the evolution 
shown in Fig.l are close to the predicted scaling behavior 
ex r^ for the trajectories in the normalized corotating 
frame. Thus we can expect a fairly moderate evolution 
for the perturbation around the L4 and L5 points dur- 
ing the inspiral of the circular parent binary, before its 
separation becomes small and the unstable modes are 
generated due to the PN effect as we see next. 

For various mass ratio /i, we followed a test particle 
initially with ga, = (72 = at r^ = 100, and measured 
the binary separation r\yu where the distance between the 
test particle and the L5 point becomes larger than 0.7 (a 
somewhat arbitrary value) in the normalized corotating 
frame. In Fig. 2, our numerical results are shown by black 
points. In addition to these results, we examined the un- 
stable separation Thu starting from finite (but small) q^ 
and qz^ and obtained similar results (see e.g. Figs. 3 and 
4) . We also studied the evolution by Newtonian equation 
of motion, and found that the binary can resonantly hold 
the test particle typically down to our termination limit 
at Tb — 6. In Fig. 2, the solid curve represents the IPN 
prediction for the binary separation rfc„ where the L4 and 
L5 points become unstable (corresponding to the degen- 
eracy condition uje = wl). At the regime \i > 0.025, the 
unstable radius Thu shows a reasonable agreement with 
the analytical result. 

Interestingly, our numerical results in Fig. 2 show a 



As we have studied so far, with RMI, a third body 
could stay around a parent BH binary deeply into the 
relativistic regime. Even though careful attention should 
be paid for interpreting our IPN results, they would pro- 
vide us with qualitative insights about the potential final 
fate of the third body. This issue would be particularly 
interesting in relation to GW observation with the Laser 
Interferometer Space Antenna (LISA) [ij. In order to 
make concrete pictures for LISA, we assume that the to- 
tal mass of a parent BH binary is M ~ 10® M©, and the 
mass of the third body is ~ 1 — IOM0. Given the stabil- 
ity condition /i < /^i = 0.03852 for a 1:1 resonance, the 
parent BH binary might be regarded as an intermediate 
mass ratio inspiral rather than an inspiral of two compa- 
rable BHs. While abundance of BHs around ^ 10'' Af© 
is not well known at present, the actual event rate of the 
former might be higher than that of the latter. 

The signal-to-noise ratio of GW associated with the 
small third body would be much lower than that of the 
parent BH binary. However, the stronger GW signal from 
the parent BH binary would enable us to easily estimate 
the basic parameters of the parent such as the mass ra- 
tio /i or the total mass M. Then, for example, from the 
estimated masses of a potential parent, we can predict 
the epoch when the L4 and L5 points become dynami- 
cally unstable (see Fig. 2). This kind of prior information 
would considerably help us to make a careful follow-up 
data analysis searching for a weaker GW signal by a third 
body. 

Statistical analyses for the final fates of RMIs with re- 
alistic initial conditions would be useful for astrophysical 
arguments, but they are far beyond the scope of this pa- 
per. Here, we would rather make qualitative discussions 
on the expected destinies. Among our numerical samples, 
ejection of a test particle from the parent binary is a fre- 





FIG. 3: Ejection of a test particle around the binary sepa- 
ration n = 24.6 after N ^ 1.1 x 10^ cycles from rt = 200. 
The ejection velocity is ~ 0.14c. The initial conditions are 
the same as Fig.l. The larger BH is at (0.027,0,0) and the 
smaller one is at (—0.973, 0, 0). They rotate counterclockwise. 



quent final state. In Fig. 3, we show the trajectory of the 
test particle evolved from Fig.l. The I/5 point becomes 
dynamically unstable around ri, ~ 24.6 and the particle 
was soon scattered by the smaller BH at (—0.973,0,0). 
It escaped from the binary with a terminal velocity of 
^ 0.14c. which is much larger than the typical kick ve- 
locity by the anisotropic GW emission [1^. For an RMI 
event, the third body can be scattered by a relativistic 
binary, and this magnitude is not surprising. The third 
body should be a neutron star or a black hole for surviv- 
ing tidal disruption during the large angle scattering by 
the smaller BH of the parent binary. 

Though an electromagnetic wave search for such a high 
velocity object would be challenging, we might get a sig- 
nature of its ejection by careful analysis of GWs. If the 
third body is a white dwarf and tidal disruption occurs, 
we might observe a short-period electromagnetic wave 
signal before the merger of the parent BH binary. This 
might help us to identify its host galaxy and perform 
cosmological studies, e.g. constraining dark energy pa- 
rameters through the relation between the redshift and 
the luminosity distance of the binary [16| . 

From our numerical samples, we expect that a plunge 
into the larger BH would be another likely scenario. Note 
that, from a distance of 0{riiu), its angular size is much 
larger than a degree scale. For a plunge into the smaller 
one, we might detect its GW signal by ground-based de- 
tectors, depending on the mass of the BH. 

Some of our numerical samples resulted in the dynami- 
cal formation of extreme- mass-ratio-inspiral (EMRI) sys- 
tems. In Fig. 4, we present the trajectory of a test particle 
around the unstable separation rt, ~ 13. After the tran- 
sitional stage shown in Fig. 4, the test particle almost 



FIG. 4: Dynamical formation of an EMRI system around Vb = 
13.0 with n — 0.0025. The initial positions are q^ — 0.004 and 
Qz = 0.08 at n = 77 (9.9 x 10* cycles to n = 13.0). 



decouples from the evolution of the parent BH binary. 
When the binary reach the innermost stable circular or- 
bit Tf, = 6, we have an eccentric EMRI system with peri- 
center distance ~ 13 and the apocenter distance ~ 50. 
GW from an EMRI is a very important target of LISA for 
directly probing highly distorted geometry around BHs, 
although we need significant effort to detect it due to 
its complicated waveform and the limitation of available 
computational power [17|. As we commented earlier, the 
basic parameters of the parent BH would be estimated 
by its stronger but simpler GW signal. Thus the subse- 
quent EMRI signal would be detected more easily than 
a blind EMRI search as usually assumed. Furthermore, 
we might, in principle, use the third body to probe the 
dynamical gravitational field caused by the merger of the 
parent BH binary to precisely measure the basic param- 
eters (e.g. mass and spin) of the merged BH. 



V. DISCUSSIONS 

Our study for RMIs is based on the IPN restricted 
three-body analysis for a circular parent binary. This 
simple treatment can be extended in diverse ways. One 
of the principle directions is to include higher-order PN 
effects and spins of BHs. A complementary approach 
would be to inject a test particle into a numerically 
evolved BH binary (for recent breakthrough see [18|). 
From Newtonian analyses we expect that modest eccen- 
tricity (e ^ 0.1) of the parent binary would not largely 
change our basic results [13|, but this should be checked 
specifically. Meanwhile, we have regarded the third body 
as a massless pointlike particle. Actually, the assurnp- 
tion for the mass would be largely relaxed for RMIs §]. 
For example, Saturn has stable co-orbital satellites (in 



a 1:1 resonance); Janus and Epimetheus whose masses 
are comparable (1:0.25) |l|. Furthermore, the internal 
structure of the third body might cause interesting as- 
trophysical effects (e.g. for a main-sequence star or even 
a gas cloud). Beyond a three-body problem, a parent 
BH binary might keep many small objects around its Li 
and L5 points, similar to the Trojan asteroids of the Sun- 
Jupiter system. In the case of a 1:1 resonance, we have 
a severe upper limit fx < fii for the restricted three-body 
problem, but a larger mass ratio might be allowed for 
resonances other than 1:1. 

In addition to these refinements for the evolution of 
RMIs, follow-up studies in related fields are worth ex- 



ploring. More detailed analyses for detectability of GW 
signatures at various stages of RMIs would be fruitful, es- 
pecially for space interferometers such as LISA. Finally, 
discussions on the formation mechanism of RMIs would 
be an attractive topic on astrophysics (see e.g. [20,] and 
references there in for related studies). 
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Appendix A: equations of motion for particle 
systems 

In this Appendix, we briefly explain the basic equa- 
tions applied in our numerical calculations for dynamical 
evolution of RMI. For the conservative IPN equations 
of motion, we use the Hamiltonian form of the Einstein- 
Infeld- Hoffman Lagrangian given for a point-particle sys - 
tem with masses nia (a: the label for particles) |9l. [2l|. 
We follow the positions Xai [i — 1,2,3: the label for 
spatial directions) in the barycentric non-rotating frame 
(xi,X2,a;3). Rather than adopting the standard conju- 
gate momentums pai, we introduce the new variables 



Pai 

rria 



(Al) 



that are convenient for dealing with restricted three-body 
problems. Our IPN Hamiltonian is expanded as 

H = Hn+ HiPN (A2) 

with the explicit forms of the Newtonian and IPN terms 

Hn = ^Z^^ias^-- 2^ — (A3) 



2 ^ rab 



Hi 



PN 



-lE^^(^l)' + l E 



marribmc 



a^b^a,c^a 



TabTa 



1 Y^ malTlb „ 2 , 7 



4 ^ Tab 



+ [riab ■ Sa){nab- Sfa)}, 



(A4) 



where Vab = \xa ~ Xb\ and Uab = [xa ~ Xb)/rab- 

For our new set of variables {xai,Sai), the canonical 
equations are modified as follows; 






1 dH dxai 1 dH 



ma dXai ' dt nia dsa 



(A5) 



These are well behaved in the limit nia — >■ for a re- 
stricted problem. Note that with the IPN terms the 
simple identity Sai = ^|f^(= Vai) does not hold, and we 
have the following correspondence at the IPN order |21| : 



-Jai ^ "at 



l"E 



b^a 



rub 

Tab 



-TtV" — {Ivu + {nab ■ Vb)nabi} ■ (A6) 

2 f-" Tab 
b^a 

In this paper, we study circular restricted three-body 
problems in which the massless third body is irrelevant 
for the evolution of the parent binary (with masses mi 
and 7712 ) • The dynamics of the circular parent binary is 
characterized by the orbital frequency wt as a function 
of the binary separation r^ = \xi — X2\, and it is pertur- 
batively expressed as 



^B — ^BN + ^BIPN 

with the Newtonian result 

1/2 



^BN 



mi + 7772 



(A7) 



(A8) 



Applying the above Hamiltonian for two particles, we 
obtain the IPN correction term as (see e.g. Q) 



^BIPN = l^BN- 



7771 + TO2 



77717772 



2rb [ (7771 + 1^2)^ 



(A9) 



Thus the motion of the parent binary can be handled 
analytically, and we can numerically follow the position 
of the third body with the modified canonical equations, 
plugging in the analytical information of the parent bi- 
nary. 

Up to now, our IPN system is conservative and this 
will provide us with a good opportunity to check the 



accuracy of our numerical code (based on a fifth order 
Runge-Kutta integration scheme [22]) for the motions of 
third bodies around the Lagrange points. For relevant 
sets of parameters (r^, /i, q^), we examined the IPN ver- 
sion of Jacobi energy derived for motions of the third 
body within the orbital plane of the parent binary [23J . 
We found that it is typically conserved by ^ 10~^ level 
for ~ 10^ orbital rotation cycles. 

Finally we summarize the dissipative effects caused by 
gravitational radiation. The first equation of (jA5|) is now 
modified as follows 



dsai _ 1 dH f dsg 
dt nia dxai \ dt 



^ > (AlO) 



where the second term is due to the radiation, and at the 
lowest order, it is given by 



l(V) 



dSm 

dt 



25 y "■^ 



(All) 



Here Z?, ■ is the fifth time derivative of the quadrupole 
moment Dij defined by (9| 



Dy = y^77ia(3xQ 



xld. 



(A12) 



For the restricted three body-problem, the moment Dij is 
determined only by the parent circular binary (a = 1,2), 
and its fifth time derivative is evaluated analytically with 
a standard adiabatic treatment. After some algebra, we 
obtain the time derivative of the separation rb of the 
parent binary as in Eq.(4), and its orbital position is 
determined straightforwardly. Then, with a method sim- 
ilar to the previous conservative case, we can numerically 
follow the massless third body. Our numerical results in 
Sec.III-IV are obtained in this manner. 



